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ABSTRACT 

Context. HD 327083 is a supergiant B[e] star that forms a binary system with an orbital semi-major axis of approximately 1.7 AU. 
Aims. Our previous observations using the VLTI and AMBER in the medium resolution K-band mode spatially resolved the envi- 
ronment of HD 327083. The continuum visibilities obtained indicate the presence of a circumbinary disc. CO bandhead emission was 
also observed. However, due to the limited spectral resolution of the previous observations, the kinematic structure of the emitting 
material could not be constrained. In this paper, we address this and probe the source of the CO emission with high spectral resolution 
and spatial precision. 

Methods. To determine the properties and kinematics of its CO emitting region, we have observed HD 327083 with high spectral 
resolution (25 & 6 km s -1 ) using AMBER and CRIRES. The observations are compared to kinematical models to constrain the source 
of the emission. 

Results. The multi-epoch AMBER spectra obtained over 5 months contain no evidence that the CO 1 st overtone emission of 
HD 327083 is variable. This indicates that the structure of the emitting region is not strongly dependent on orbital phase. It is shown 
that the CO bandhead emission can be reproduced using a model of a Keplerian disc with an inclination and size consistent with our 
previous VLTI observations. The model is compared to AMBER differential phase measurements, which have a precision as high 
as ~30 yuas. A differential phase signal corresponding to 0.15 mas (~5 cr) is seen over the bandhead emission, which is in excellent 
agreement with the model that fits the CRIRES observations. In comparison, a model of an equatorial outflow, as envisaged in the 
standard sgB[e] scenario, does not reproduce the observations well. 

Conclusions. We present a direct test of the circumstellar kinematics of the binary sgB[e] star HD 327083 using both spatial and 
spectral information. The excellent agreement between the disc model and observations in the spatial and spectral domains is com- 
pelling evidence that the CO bandhead emission of HD 327083 originates in a circumbinary Keplerian disc. In contrast, the model 
of an equatorial outflow cannot reproduce the observations well. This suggests that the standard sgB[e] scenario is not applicable to 
HD 327083, which supports the hypothesis that the B[e] behaviour of HD 327083 is due to binarity. 

Key words. Stars: circumstellar material - Stars: early type - Stars: emission-line, Be - Stars: mass-loss- Stars: individual: HD 
327083 

1. Introduction as gamma ray bursts may be constructed in the sgB[e] stage. 

Consequently, it is important to study the circumstellar geome- 
Massive stars play a crucial role in multiple areas of astro- of sgfi[e] starfj and determine what is reS p 0n sible for their 
physics. As a result, a complete understanding of many astro- mass j oss 
physical phenomena requires an understanding of the evolution 
H ' of massive stars. In turn, this requires knowle dge of how mas- With this in mind, we recently studied the sgB[e] star 
■ ■' sive stars lose mass (see e.g. iPuls et al.ll2008l Supergiant B[e] HD 327083 with the VLT I and AMBER in the med ium reso- 
(sgB[e]) stars are important objects in this regard as they are lution ^-band mode (see |Wheelwright et alj|2012bj hereafter 
massive objects in a late evolutiona ry stage which also exhibit W12). High spatial resolution observations were required to 
signs of enhanced mass loss (see e.g. lLamers eUdll998l In ad- P robe the geometry of the object's circumstellar environment on 
dition, it has been suggested that there could be an evolution- milli-arcsecond scales. In turn, this was necessary to elucidate 
ary lin k between sgB [ el stars and the luminous blue variables the "bracket" behaviour of HD 327083, for which there have 
(LBVs, IZickgrafll20Q6l) , In turn, LB Vs are generally thought to been several hypotheses. On the basis of a comparison between 
be the precursors of the Wolf -Rayet stars ; the final phase in mas- optical sp ectroscopy and a NLTE mod el of an expanding at- 
sive star evolution (see e.g. IVinkll2009l) . Therefore, the asym- mosphere,|Machado & de Araujo| (|2003|) suggest that the object 
metric geometries that are required to produce phenomena such ma Y be close to the Lfi V phase. This would make it a key object 

to study the link, if any, between th e sgB[e] and LBVs. However, 

* Based on data obtained at the European Southern Observatory, there is an alternative hypothesis. |Miroshnichenko et al.| (|2003|, 
Paranal, Chile via the observing programmes 385.D-0513 and 087.D- hereafter M2003) detected an unresolved binary companion via 
0889. radial velocity variations. Based on the relatively short period of 
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the variations, ~6 months, these authors suggest that the system 
is close enough to interact. In this case, the material traced by 
the object's infrared excess may be mass lost as the result of bi- 
nary interactions, as proposed fo r othe r objects showing the BT el 
phenomenon (see lMillour et al.ll2009l I201U iKraus et al.1120121) . 

We now briefly summarise the results from W12. Our obser- 
vations with the VLTI and AMBER spatially resolved the ma- 
terial responsible for the /f-band excess of HD 327083. Using 
simple geometrical models, we found that the continuum visibil- 
ities could be reproduced by an elongated ring with a central ra- 
dius of approximately 6.6 AU. The location of this ring is consis- 
tent with the expected dust sublimation radius of a star with the 
stellar parameters reported by M2003; 6 < R su b < 25 AU u sing 
the standard equation (see iMonnier & Millan-Gabetll2002l) . A 
non-zero closure phase was observed over the continuum emis- 
sion and this was attributed to the binary companion. The closure 
phase suggested that the companion should be located within the 
elongated ring. Therefore, the geometrical model was associated 
with the inner rim of a dusty, circumbinary disc. Such a config- 
uration is more rem iniscent of mass los t as the result of binary 
interactions (see e.g. lMillour et al.l201 ll) than intrinsic mass loss 
via an equator ial outflow, as depicte d in the standard sgB[e] sce- 
nario (see e.g. lZickgraf et al.lll985l) . Therefore, the observations 
to date favour the interacting binary hypothesis. 

Probing the mass loss process of the system further requires 
constraining the kinematics of the circumbinary material and 
the conditions that give rise to the occurrence of the B[e] phe- 
nomenon. Such constraints can be determined from spectrally 
resolved observations of CO bandhead emission. Therefore, we 
obtained new high spectral resolution (25 & 6 kms -1 ) obser- 
vations of HD 327083's CO bandhead emission with AMBER 
and CRIRES. Here we present these observations, which provide 
sub-milli-arcsecond precision constraints on the geometry of the 
emitting material. Th e CO bandhead e missio n of HD 327083, 
initially observed by iMcGregor et al.l (119881) . is spectrally re- 
solved for the first time and we observe a differential phase sig- 
nature over the emission. This paper presents the observations 
and modelling of the CO emitting material, and it is structured as 
follows. In Sect. [2] we describe the observations and subsequent 
data reduction. We present the observational results in Sect. [3] 
These are then analysed and modelled in Sect. [4] and the impli- 
cations are discussed in Sect. The conclusions are presented 
in Sect. [6] 

2. Observations and data reduction 

2.1. CRIRES 

HD 327083 was observed with the near-IR cryogenic high res- 
olution spectrograph CRIRES of the VLT (iKaeufl et alj|2004 
during the night of June 28, 2010. The target was included in 
our survey of sixteen Galactic, A'-band bright sgB[e] candidates 
and unclassified B[e] stars. The aim of the project is to study in 
detail the kinematic signature of the CO first overtone bandhead 
of B[e] stars in general (Muratore et al. in prep.). Here we make 
use of the data for this object. Observations were executed with 
a slit-width of 0.4" delivering a spectral resolution of approx- 
imately 50,000 (or 6kms~ 1 ). The data were obtained applying 
a standard nodding on-slit strategy to remove sky and detector 
glow, without AO correction and applying a small jittering offset 
at each nod position to better remove bad pixels. The DIT was 
chosen bearing in mind the non-linearity regime of the CRIRES 
InSb Aladdin III detectors. A telluric standard was taken imme- 
diately after the science targets, close in airmass. The seeing con- 
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Fig. 1. The projected baselines of the AMBER observations of 
HD 327083. 

ditions were poor as the seeing varied between 1.75 and 2.25". 
However, the brightness of the target, K =3.3 mag., ensured the 
data are still of high quality. Data reduction was conducted in a 
standard fashion utilising the ESO provided CRIRES pipeline. 

2.2. AMBER 

HD 327083 was observed with the VLTI and the 3-beam com- 
biner instrument AMBER dPetrov et al. 2007|) in the high spec- 
tral resolution /if-band mode centred on 2.288 fim. This setting 
delivers a spectral resolution of R =12,000 and a wavelength 
range of 2.265 to 2.31 1 fim. FINITO was used to provide fringe 
tracking. Observations were conducted on 6 occasions between 
May 2 nd 201 1 and September 22 nd 201 1 . The observations were 
carried out employing the 1.8m Auxiliary Telescopes (ATs). A 
log of the observations is presented in Table [T]and the projected 
baselines are presented in Fig.Q] In general, short baselines were 
used as our previous observations demonstrate that this object is 
relatively extended. The seeing conditions were typically good 
(the mean seeing was ~0.75"), which was important to facili- 
tate these high spectral resolution observations. Observations of 
HD 327083 took place between observations of calibrator ob- 
jects (HD 151078 and HD 161068). 

The AMBER observations were reduced in a st andard fash- 
ion using the amdlib s oftware version 3.0.3 (see iTatulli et"a"D 
120071: IChelli e7ai]l2009h . The processes involved, i.e. determin- 
ing the pixel-to-visibility-matrix and converting the observed 
fringe patterns to measurements of the coherent flux, are stan- 
dard operations with this software and are not described here. 
Once the interferometric observables are calculated, the general 
approach is to subject each block of frames to certain selec- 
tion criteria before averaging the selected frames into a single 
measurement. Since we employed the high spectral resolution 
mode of AMBER with a long integration, and the data exhibit a 
range of frame signal-to-noise-ratios, it is challenging to accu- 
rately calibrate the visibilities. Therefore, we do not make use 
of absolute visibilities in this paper. As we are not concerned 
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Table 1. The log of the AMBER observations. 



Date 


Stations 


Seeing 


Coherence time 


1J1 1 


Baseline 


PA 






c) 


(ms) 


(s) 


(m) 


o 


02-05-2011 


H0-I1-G1 


0.46 


5.3 


12.0 


40.7/38.6/65.7 


169.9/58.1/23.0 


07-05-2011 


D0-A1-C1 


0.68 


4.1 


12.0 


34.7/15.7/21.8 


45.4/71.7/26.9 


12-06-2011 


D0-B2-C1 


0.88 


4.6 


12.0 


33.2/11.0/22.1 


22.2/22.1/22.2 


27-07-2011 


K0-I1-G1 


1.03 


2.1 


12.0 


44.2/43.6/84.9 


19.3/49.0/34.0 


04-08-2011 


K0-I1-G1 


0.74 


2.1 


12.0 


45.1/45.9/88.1 


11.1/39.5/25.4 


22-09-2011 


K0-I1-G1 


0.77 


3.1 


12.0 


38.2/29.5/65.1 


33.9/66.2/47.9 



Notes. The baseline lengths and position angles (PAs) represent time-averaged values for a given night. 



with temporal variations in the transfer function, we merge all 
consecutive observations before performing frame selection and 
averaging. This provides a higher resultant signal-to-noise-ratio 
(SNR) (see e.g. lStefl et aT1l2009h . A high frame selection rate, 80 
per cent, was used as this generally provides the highest quality 
phase information and our analysis in this paper concentrates on 
these measurements. 

Wavelength calibration of the AMBER data was conducted 
using the telluric lines present in the spectra. To provide a ref- 
erence wavelength for the telluric lines observed, we utilised a 
high resolution (R =40,000) spectrum of the Earth's telluric fea- 
tures made available by the NSO/Kitt Peak ObservatoryQ. The 
position of the telluric lines was fit with a 2 nd order polynomial 
which was used to provide the wavelength calibration. The mean 
rms of the solutions was found to be approximately 4 kms -1 , 
which was deemed acceptable. Only the spectra of HD 327083 
were used to provide wavelength calibration as the spectra of the 
standards exhibit CO absorption features, which made identifi- 
cation of the telluric lines less secure. 

On several occasions, the closure phases of the calibrator 
objects were non-zero. In an attempt to correct this, we sub- 
tracted the mean of the calibrator closure phases from those of 
HD 327083. Since we merged the OI files prior to frame selec- 
tion and averaging, the closure phase data calibration could be 
compromised by temporal variations in the instrumental closure 
phase. In general, the closure phases of the two calibrators are 
consistent, indicating this is not an issue. However, the calibra- 
tor closure phases for the observations conducted in September 
were inconsistent by approximately 5°. This translates to a sys- 
tematic uncertainty in the final calibrated measurement. 

The two calibrator stars, HD 151078 and HD 161068, were 
also used as telluric standards. For each night of observing, we 
attempted to remove the telluric absorption features in the spec- 
trum of HD 327083. Before this could be done, the CO ab- 
sorption features in t he spectra of the standard had to be re- 
moved (see e.g. W12. lTatulli et al.ll2008l) . This was done by di- 
viding the standard spectra with a telluric corrected spectrum 
of a template star with the same spectral type. High-resolution 
(R~ 18,000), /T-band template spectra were obtained from the 
library of 25 standard stars observed with GNIRS provided by 
GEMINI observatoiyS For HD 161058 (K4/K5II0) the spec- 
trum of HD 78541, also of spectral type K4/5III was used. In the 
case of HD 151078 (K0III) the average spectrum of HD 74137 
and HD 28 (both K0III) was used. 



1 http://www.eso.org/sci/facilities/paranal/instruments/isaac/tools/ 
spectroscopic_standards.html#Telluric 

2 http://www.gemini.edu/sciops/instruments/nearir-resources/ 
spectral-templates/rl 8000/ 

3 from SIMBAD: http://simbad.u-strasbg.fr/simbad/ 



3. Observational results 

The AMBER and CRIRES observations spectrally resolve the 
CO bandhead emission of HD 327083 for the first time (the 
spectra are presented in Sects 14. Il and l4.2l respectivelv). The high 
spectral resolution of CRIRES (~6 km s~') and the high SNR of 
the data (~200) provide a unique view of the bandhead spectral 
profile. The bandhead clearly displays a blue wing, the so-called 
"blue shoulder", which can indicate that the emission originates 
in a r otating medium such as a disc (see e.g. lWheelwright et alj 
1201 Oi hereafter W10, and references therein). No clear signa- 
ture is seen over the bandhead emission in the differential visi- 
bilities. However, we report the first detection of a signal across 
CO bandhead emission in AMBER differential phase observa- 
tions (this is presented in Sect. 14.3b . This allows us to probe the 
emitting region with a precision as high as 30 /^mas. 

Turning to the closure phase measurements (presented in 
Sect. l4.4t , no clear signature is seen over the CO bandhead emis- 
sion. In general, the average closure phases are close to zero; 
which suggests that the environment of HD 327083 is relatively 
symmetric. However, the observations in July suggest a shift to 
non-zero values, which indicates an increase in the asymmetry 
of the environment of HD 327083. In the following sections we 
analyse these points in more detail to constrain the source of the 
CO emission and the nature of the mass loss of HD 327083. 

4. Analysis and modelling 

4.1. CO bandhead spectra: assessment of variability 

Our AMBER observations of HD 327083 span a period of ap- 
proximately 5 months, which is more than 80 per cent of the 
orbital period estimated by M2003. Therefore, we can use the 
observed flux spectra to assess whether the CO emission is vari- 
able, which might indicate mass transfer at particular orbital 
phases. To determine whether there are significant differences 
between the average spectra from each night o f observation, w e 
use the temporal variance spectrum (TVS, see lFullertonlll990l) . 
The TVS can be expressed by: 

where N represents the number of spectra, Fj(A) denotes the 
individual spectra, F av (A) is the average spectrum and <x, is the 
uncertainty at each point of the individual spectra. Here, we used 
the rms in the continuum as a constant uncertainty across the 
individual spectra. 

Following the example of lOudmaiier & Bakkerl (1 19941) . we 
calculated the temporal sigma spectrum, TSS, which is given by 
TSS = yTVS. This value is expressed in terms of the noise 



3 



H.E. Wheelwright et al,; Tracing the origin of the CO bandhead emission of HD 327083 



jl.2 
jl.O 



£l.O 
$0.8 



- s% 

^~ J v Average 1 


yuy Ti Ai ysj y ^ 

elluric corrected- 




Uncorrected- 









2.290 



2.295 

X (fJ-TTl) 



2.300 



Fig. 2. The upper panel presents the time-averaged, telluric cor- 
rected AMBER spectrum. The central panel contains the un- 
processed (not corrected for telluric absorption) AMBER spec- 
trum. The shaded region represents the scatter of the six aver- 
age spectra constructed for each observing date. The lower panel 
presents the temporal sigma spectrum. 



level. Therefore, if there are no significant deviations between 
the individual spectra, this value will be approximately unity. On 
the other hand, significant variations will be associated with a 
large TSS value. The average AMBER spectrum of HD 327083 
(uncorrected for telluric absorption) and the associated TSS is 
shown in Fig. [2] HD 327083 was observed at different air-masses 
each of the six times it was observed. As a result, the strength of 
the telluric absorption features in each spectra is different, al- 
though, as shown in Fig. [2] the differences are relatively small. 
The telluric lines in the average spectrum are associated with 
peaks in the TSS. This confirms that this approach is sensitive to 
variability. The bandhead itself, however, exhibits no such fea- 
ture. Therefore, we conclude that these data contain no evidence 
of variability in the CO bandhead emission. 



We note that the final telluric corrected AMBER and 
CRIRES spectra are essentially consistent, which is to be ex- 
pected if the emission is not variable. However, comparison of 
the two suggests that the telluric correction in the AMBER spec- 
trum is not perfect. This is likely because the standards were 
chosen as interferometric calibrators rather than telluric stan- 
dards and thus they do not necessarily provide the best air-mass 
match. From here-on, as far as flux spectra are concerned, we 
concentrate on the high spectral resolution CRIRES spectrum. 

4.2. Modelling the CO flux spectrum 

With the aim of constraining the origin and kinematics of the 
bandhead emission, the observed spectral profile is fit with a 
kinematical model. Specifically, to evaluate the disc scenario, we 
attempt to fit the CO bandhead with a model of a rotating disc. 
We have already used a similar model to fit the CO emission of 
massive young stellar objects in other publications (W10, Ilee et 
al. in prep.), and we thus only describe the methodology briefly 
here. 

A simple, geometrically flat disc in Keplerian rotation is con- 
structed. The excitation temperature and surface number density 
decrease with increasing radius according to power laws. The 
temperature decreases with R~° J5 and the surface number den- 
sity falls off with R l 5 , in line with standard disc theory (see 



e.g. lKraus et al. 2000). We note that the temperature is the exci- 
tation temperature used to determine the relative populations of 
the CO ro-vibrational energy levels. Therefore, it may not neces- 
sarily follow the thermal temperature distribution. F or example, 
it is possible the vi brational levels are not in LTE (iNaiita et al.l 
1 19961 : iMartinl 1 1 997h . The excitation temperature distribution is 
difficult to constrain with only one bandhead p rofile as it can 
vary significantly (see e.g. iBerthoud et al1l2007l) . Therefore, we 
fix the temperature distribution as described. 

The temperature and surface density at the inner radius, R^, 
are free parameters. The CO emission is assumed to originate 
in the region of the disc where 1000 K < T < 5000 K. In the 
case of HD 327083, we set the inclination of the disc to be 48.5°, 
as determined by W12. In most cases, it would be assumed that 
the disc extends from the stellar surface. However, in the case 
of HD 327083, it is possible that the disc is circumbinary and 
the inner disc is cleared. Therefore, we allow the inner rim of 
the disc to vary. The free parameters of this particular model are 
thus: the inner rim of the disc (technically the radial distance 
of the smallest annulus), the temperature at R m , the intrinsic 
line-width and the surface number density at /?;„. The disc was 
split into radial and azimuthal cells and the CO emis sion of each 
cell w as calculated according to the methodology of lKraus et al.l 
(120001) . Finally, once the individual cell spectra were calculated, 
they were summed to create the total spectrum. This was then 
smoothed to match the spectral resolution of the observations. 

This relatively simple approach is standard in such studies as 
there are no sophisticated models including CO bandhead emis- 
sion available. As the code is analytic, it is readily used in a 
fitting scheme. The best-fitting model was determined using the 
downhill simplex algorithm (implemented in idl as the amoeba 
routine). Once the best fit model was found, monochromatic im- 
ages were created, which were subsequently compared to the 
interferometric data. The uncertainties associated with the final 
parameters were determined in the following way. A 2D reduced 
X 2 map was constructed for each set of parameters, centred on 
the best fitting model. The uncertainties in the best fitting pa- 
rameters were estimated by determining the maximum parame- 
ter change that results in a chang e in x 2 of 1 . As discussed in 
iDomiciano de Souza et al.l (1201 ll) . this likely results in an over- 
estimation of the uncertainty (provided the data points are un- 
correlated). To investigate the possibility of additional fits with 
different inner radii, we constructed 2D x 1 maps centred on the 
best fitting model but displaced in the direction of the inner ra- 
dius by -5 and +5 R* (we concentrated on the inner radius as 
this influences the interpretation of the data, as discussed in Sect. 
[5). In both cases, the minimum^ 2 was found to be greater than 
that of the model presented in Table [2] Therefore, we conclude 
that the best fitting model does indeed correspond to the global 
X 2 minimum. 

The model spectrum that is found to best reproduce the CO 
bandhead profile is shown in the top panel of Fig. [3] The param- 
eters of the best fitting model and the adopted stellar parameters 
are presented in Table [2] In general, the model provides an ex- 
cellent match to the features of the bandhead. Indeed, the fit is 
remarkably good considering the limited number of free param- 
eters and the simplicity of the model. Therefore, we conclude 
that a model of a Keplerian disc with an inner radius of approxi- 
mately 3 AU and a relatively low initial temperature, 7Ao~ 1700 K, 
appears to be a fitting representation of the CO emitting region 
associated with HD 327083. 

However, the possibility remains that the source of the CO 
emission is not rotationally dominated. Indeed, the favoured 
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model of the circumstellar environment of sgB[e] stars features 
an equatorial outflow. To investigate whether this scenario is 
applicable to HD 327083, we perform an additional fit to the 
CRIRES spectrum with an outflow model. 

The standard model for the sgB[e] environment is a dense 
equatorial outflow and a faster wind concentrated in the polar 
regions. The equatorial outflow is thought to be the result of 
rapid rotation and the bi-stability jump leading t o enhanced line 
driven mass loss in the equatorial region (see e.g. lPelupessv et ail 
In modelling the emission fr om an equatori al outflow, 
we adopt the analytical treatment of iPorterl (120031) who com- 
pared the SEDs of disc and outflow models to observations of 
the sgB[e] star R126. Specifically, the CO emission is modelled 
as originating in a flat structure containing both rotation and out- 
flow. The velocities are given by: 

v Rot = yjGM*/R+\j^-\ (2) 
and 

VOut = VSound + (Voc ~ V Sound ) 1 1 - -y J (3) 

which represent angular momentum conserving rotation and the 
velocity structure in a line driven wind respectively. 

We fit the CRIRES spectrum with a geometrically flat, equa- 
torial outflow model featuring both these velocity laws by treat- 
ing the terminal velocity, line-width, initial number density and 
temperature as free parameters and setting = 1 R+ (since the 
wind is launched from the star) and fi = 0.8 (which is appro- 
priate for a hot star wind). The inclination was also set, as in 
the disc modelling. The temperature and number density distri- 
butions follow the same power laws as in the disc model, which 
is appropriate for this analytic treatment (see e.g. t he discussio n 
on the heating of circumstellar discs and winds in Porter 2003). 
As before, the best fitting model was found using the downhill 
simplex algorithm. 

The best fitting spectral profile of the equatorial outflow 
model is also shown in the upper panel of Fig. [3] The best fit- 
ting parameters are: N = 2.7 x 10 22 crrT 2 , T = 19624 K, 
Linewidth = 6.8 kms" 1 and V m = 69.4 kms" 1 . As can be seen, 
the outflow model can approximately reproduce the CRIRES 
spectrum, although the quality of the fit is noticeably worse than 
that provided by the disc model. This is reflected in the large 
reduced^ 2 , which is ~100. Therefore, the^ 2 values favour the 
Keplerian model. Nonetheless, since the spectrum is approxi- 
mately recreated, it is not possible to confidently discard the out- 
flow model based on the CO bandhead profile alone. However, 
the addition of our AMBER data provides another opportunity 
to differentiate between the two models. This is described in the 
following section. 



4.3. Differential visibilities and phases 

The spectral modelling of the previous section provides detailed 
predictions regarding the spatial distribution of the CO emis- 
sion. Therefore, the AMBER data offer a unique test of the 
spectral models that fit the CO bandhead. No clear signature is 
seen in the differential visibilities across the bandhead emission. 
Considering the uncertainties in the visibilities (typically 10 per 
cent) and the low line-to-continuum flux ratio (~0.2), this is con- 
sistent with the CO bandhead emission originating close to, and 



Table 2. The parameters of the CO emitting Keplerian disc 
model that best fits the CRIRES spectrum. 



Parameter 


Vtiluc 


INOLcS 


Cent. Mass 


25 M Q 


Set 1 




25 R 


Set 2 


i 


45.8 


Set 3 


N 


0.4 ± 1.8 x 10 24 cm" 2 






25.5+2.3^ 


~3 AU, 2 mas 


Linewidth 


3.2+3.2 km s" 1 




To 


1721±540 K 




x 2 


3.96 





Notes. 1: Based on an early B type primary and an F type companion 
(M2003), 2: This sets the minimum possible inner disc radius, based on 
the values reported by M2003, 3: Set by our previous AMBER obser- 
vations (W12) and assuming that the CO emission is co-planar with the 
continuum emission. The outer radius is simply the radius at which the 
temperature falls below 1000 K, and in this case is approximately 6 AU. 



perhaps interior to, the continuum emission (see the discussion 
in W12). Our spectral models place the CO material interior to, 
but extending to, the continuum emitting region. Therefore, the 
best fitting models are consistent with the observed visibilities 

A second, more stringent test is provided by the AMBER 
differential phase measurements. Differential phases are sensi- 
tive to photo-centre displacements associated with asymmetri c 
sources of flux, such as rotating discs (see e.g. lLachaumel200"3h . 
Therefore, we compare the differential phase measurements to 
the predictions of the models that fit the CRIRES data. To cal- 
culate the photo-centre offsets of the models, it is assumed that 
the continuum is centro-symmetric. Most of the closure phases 
measured are close to zero (see Sect. 14.4) . indicating that the en- 
vironment of HD 327083 is relatively symmetric, justifying this 
approach. 

To improve the SNR of the data, we re-bin the differential 
phases by a factor of 3 (approximately a resolution element). 
The minimum uncertainty was ~30 //as. This uncertainty corre- 
sponds to the data obtained with the 34.7 m D0-A1 baseline used 
on 7 th May 201 1 . These data alone are precise enough to test the 
models at a significant level. The errors on the other baselines 
are generally a factor of two or more greater, so we concentrate 
on the data obtained with this baseline only. 

A comparison between the predicted photo-centre offsets 
and that calculated from the differential phase measurements us- 
ing the aforementioned baseline is presented in the lower panel 
of Fig. [3] The offsets calculated from the differential phases 
across the other baselines are presented in Fig. 2] The combi- 
nation of high spatial and spectral precision provides stringent 
constraints on the kinematics of emission line regions. As can 
be seen, the predicted offsets of the best fitting disc and outflow 
models are different and thus we can use the spatial information 
provided by the differential phases to differentiate between the 
competing models. 

We focus on the best fitting disc model first. As can be seen 
in Fig. [3] the maximum predicted offset occurs over the band- 
head shoulder where many individual ro-vibrational lines are 
super-imposed (as discussed in W10). The offsets over the indi- 
vidual ro-vibrational lines are clearly too small (<0. 1 mas) to be 
detected here. However, the predicted offset over the bandhead 
shoulder, ~0.15 mas, is large enough to detect with the current 
data, which has an noise level of ~30 //as. The observations dis- 
play a displacement of approximately 0.15 mas over the band- 
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head shoulder, in excellent agreement with the prediction of the 
model that fits the CRIRES spectrum. The significance of the 
peak pixel of the observed signature is approximately 5 cr. The 
agreement between the amplitude of the model and observed sig- 
natures was optimised by slightly changing the position angle of 
the disc (by 10° so that the difference between the position an- 
gle of the major axis of the disc and the baseline was ~42°). This 
slight difference in PA could represent uncertainty in the previ- 
ous estimate. However, it should be noted that the distance to HD 
327083 has an uncertainty of approximately 30 per cent, which 
translates to a similar uncertainty in the model photocentre shift. 
Therefore, the agreement between the model and observations 
could also be optimised by slightly reducing the distance to HD 
327083. As a result, the minor change in PA is not considered 
significant. 

We emphasise that bar a modest adjustment of the PA, the 
observed differential phase signature was not fit. The model sig- 
nature was predicted by the fit to the CRIRES spectrum. The 
excellent agreement between the model and the data in both the 
spectral and spatial domains strongly suggests the model is an 
accurate representation of the CO emitting region of HD 327083. 
We note that if the profile were fit with rotation that decreased 
more rapidly with radius than the Keplerian case, for example 
angular momentum conserving rotation in which v oc r _1 , the 
resultant positional offse t would be smaller than th at observed 
and vice versa (see e.g. IWheelwright et alj|2012al) . Therefore, 
the agreement between the model and data is compelling evi- 
dence that the CO bandhead emission traces a Keplerian disc. 

Unlike the disc model, the outflow model cannot reproduce 
the data well. As shown, the photo-centre offset associated with 
the best fitting outflow model is smaller than that of the best fit- 
ting disc model. This is because the velocity structure and bright- 
ness distribution of the two best fitting models are very different. 
This is shown in Fig. [5] The CO emitting region of the outflow 
is clearly smaller than in the disc case. As a result, the spatial 
signature of the outflow is smaller than that of the disc and does 
not reproduce the differential phase shift seen (see Fig. |3). The 
discrepancy between the model and observations is too large to 
be accounted for by a modest change in PA or distance. We con- 
clude that the large x 1 of the spectral fit and the poor reproduc- 
tion of the photo-centre offset imply that the source of the CO 
emission is not an equatorial outflow and that these data favour 
the disc scenario. 



4.4. Closure phases 

The AMBER closure phase signa l, sampled at different phase s 
of the binary period (~6 months, Mi roshnichenko et al.ll2.Q03T) . 
has the potential to trace changes in the degree of symmetry in 
the environment of HD 327083. This might be expected due 
to changes in the binary separation and PA and would allow 
us to constrain the binary orbit to a much higher degree. The 
wavelength-averaged closure phase measurements are presented 
in Fig. [6] All bar one of the closure phase measurements indicate, 
or are consistent with, a symmetric environment. Therefore, we 
do not definitively detect the binary companion. To investigate 
the geometrical constraints supplied by the closure phase obser- 
vations, we fit the observations with two geometric models: tha t 
of a binary and the skewed ring model of iMonnier et al 1 (12006b . 
The binary model has only two free parameters, the binary PA 
and separation. The two stellar components are treated as un- 
resolved uniform discs and the binary flux ratio was set to 3.5, 
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Fig. 6. The wavelength-averaged calibrated closure phase for 
each observing date (top panel) and an example of the closure 
phase data (data obtained on the 27 th of June, lower panel). The 
data have been re-binned by approximately twice the resolution 
element. 



as discussed in W12. The parameters of the skewed ring model 
are: the ring radius and elongation, the ring flux, the ring PA, the 
degree of skewness and the PA of the skewed distribution. The 
properties such as the ring radius, inclination and flux were set 
to those of the ring found by W12 to represent the continuum 
emission of this object. Therefore our skewed ring model also 
has 2 unknowns: the degree of skewness and the PA of the skew. 

It was found that neither model could be excluded. 
Ultimately, the observations hint at an increase and subsequent 
decrease in the asymmetry of the environment of HD 327083. If 
this represents a change in the amount of asymmetry in the struc- 
ture resolved in W12, its degree of asymmetry must increase on 
a timescale of several months. This appears unlikely given the 
size of the structure (approximately 5 AU). In addition, we note 
that the CO bandhead emission is located close to and interior to 
the continuum emission. Therefore, if the continuum emission 
varied, perhaps as a result of intermittent mass loss or accretion, 
variability in the CO emission would be expected. Since the CO 
emission is not variable, it appears that the continuum emission 
originates in a stable structure. Consequently, we attribute our 
tentative detection of a change in closure phase to motion of the 
binary system. A conclusive confirmation of the interferometric 
binary detection may require model independent imaging with, 
for example, the PIONIER instrument (e.g. as demonstrated by 
iBlind etal.ll201ll) . 
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Fig. 3. A comparison between the observations and the two physically different best fitting models explored in the text. The upper 
panels present the average CRIRES spectrum and the profile of the best fitting model. In the lower panels, the observed photo-centre 
offsets calculated using AMBER differential phase measurements from the baseline with the lowest noise level (D0-A1 used on the 
7 th of May) are shown. For each model that fits the CRIRES spectrum we over-plot the predicted signature for the selected baseline. 
The model signatures are smoothed to the spectral resolution of AMBER and are shown re-binned using the size of the AMBER 
pixels and with a finer sampling. 



5. Discussion 

This paper presents high spectral resolution (25 & 6 kms -1 ) 
AMBER and CRIRES observations of the binary supergiant 
B[e] star HD 327083. The combination of information on AU 
sized scales provided by the AMBER data and the velocity re- 
solved information provided by the CRIRES data offers new in- 
sights into the close environment of this object. In particular, 
we show that the data can be reproduced well with a Keplerian 
disc model. On the contrary, the scenario of an equatorial out- 
flow appears unlikely since such a model cannot reproduce the 
observations. 

These observations, and the disc model that best fits the 
CRIRES spectrum, confirm several previous suggestions that 
were based on spectrally unresolved interferometric observa- 
tions. In W12 it was remarked that there was no clear visibility 
signature across the CO overtone emission. On this basis, it was 
suggested that the CO emission originates from a region with 
the same elongation, and thus inclination, as the near-IR contin- 
uum. Furthermore, in W 12 it was suggested that the CO material 
is located interior to the continuum emission. The geometry of 
the CO bandhead emission predicted by the spectral fit to the 
CRIRES spectrum is consistent with both these hypotheses. The 
near-IR continuum itself is thought to be dominated by thermal 
emission of warm dust. Therefore, the location of the CO emit- 
ting region is consistent with the presence of a gaseous disc in- 
terior to the near-IR continuum emitting region (which has an 
inner radius of approximately 5 AU). This supports the associ- 
ation of the continuum emitting region with the inner rim of a 
dusty disc. 

The most compelling support for the disc scenario is pro- 
vided by the kinematic element of the modelling. The excel- 
lent agreement between the model spectral profile and the data 
proves that the kinematics of the emitting material are entirely 
consistent with Keplerian motion. Finally, we note that the 



AMBER differential phase signature provides a third argument 
in terms of photocentric displacement on scales of 0.15 mas. 
The observed signature is found to be fully consistent with the 
predictions of the best fitting model. The combination of high 
spatial and spectral precision provides stringent constraints on 
the kinematics of line emission regions. Therefore, the excellent 
agreement between the data and disc model in both the spatial 
and spectral domains is compelling evidence that the CO emis- 
sion originates in a Keplerian disc. 

As far as the nature of this disc is concerned, we observe 
the following. We note that the relatively low inner temperature 
of the disc suggests there is an absence of hot material close to 
the central star. Consequently, this is consistent with the pres- 
ence of an inner gap (see Fig. |5}. The derived inner radius of 
the CO emitting disc is ~3 AU, which is larger than the (some- 
what uncertain) binary separation of a ~1.7 AU (M2003). In 
other words, the model presented here confirms the suggestion 
of W12 that the HD 327083 system is surrounded by a circumbi- 
nary disc. The fact that the gaseous disc does not extend closer 
to the central star, where its inner rim temperature would be 
~5000 K, provides further insight into the mass loss of this sys- 
tem. Specifically, the classification of HD 327083 as an interact- 
ing f3 Lyrae type binary (M2003) can most likely be discounted, 
because of the absence of a circumstellar accretion disc. Also, 
the sgB[e] model appears unfitting for this system (H D 327083 
is listed as a se B[e] candidate by iKraus etaT]l2009h . The CO 
bandhead emission of such objects is expected to trace an equa- 
torial outflow. Since we show that this model does not reproduce 
the data well, it appears unlikely that this object can be described 
by the standard sgB[e] picture. This supports the hypothesis that 
the B[e] behaviour of HD 327083 is related to its binarity. 

This conclusion raises the question of whether there are any 
other systems similar to HD 327083. The CO emission con- 
tains two clues regarding the origin of the material surrounding 
HD 327083. First, we reiterate that the material is confined to a 
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Fig. 4. Average differential phases converted into photo-centre shifts for each night of observing. The predictions of the disc model 
that fits the CRIRES spectrum are over-plotted. 



circumbinary disc. Secondly, we note that the gaseous circumbi- 
nary disc exhibits no evidence of variability. CO emission traces 
hot dense gas in the inner-most regions of discs and is sensi- 
tive to the surface density of the emitting region. Therefore, it is 
likely that the CO emission of HD 327083 would be sensitive to 
varying mass injection into the circumbinary disc. Consequently, 
the lack of variability suggests that the circumbinary disc is a rel- 
atively stable structure. 

Considering these points, the binary HD 327083 appears 
very si milar to the b i nary A [el star HD 62623 recently stud- 
ied bv IMillour et al.l (|2011|) . HD 62623 is also a binary and 
IMillour et al.l ( 20TT1) demonstrate that a circumbinary disc sur- 
rounds the system . To explain the characteristics of this system, 
iPlets et al] d 1995b proposed that angular momentum transfer oc- 
curs through the L2 Lagrangian point of the system. As a result, 
mass lost from the primary supergiant by its radiatively driven 
wind is deposited into a stable dens e circumbinary disc (which 
is consistent with the observations of lMillour et al.ll20"l"hV If this 
scenario applies to HD 327083, the mass injection into the disc 
would be less dependent on orbital phase than the Roche Lobe 
overflow at periastron suggested by M2003. Therefore, this sce- 
nario could explain the lack of variability in the CO emission. 
Furthermore, this scenario also explains the presence of a cir- 
cumbinary disc. Consequently, the environment of HD 327083 



is entirely consistent with this scenario. As a result, HD 327083 
could belong to the class of massive binaries which feature cir- 
cumbinary discs thought to be due to the channeling of the pri- 
mary's wind through the L2 point of the system (cu rrently exem- 
plified by HD 62623 and p ossibly also v Sgr, see IMillour et al.l 
l2^lllBonneau et al.ll20TTh . 



As a final remark, we note that this paper demonstrates the 
diagnostic potential of high spectral resolution observations of 
CO bandhead emission coupled with information on small an- 
gular scales as provided by VLTI observations. This approach 
offers a means to directly probe the circumstellar kinematics 
of sgB[e] stars. Recent observations have found additional ev- 
idence for Keplerian structures around other sgB[e] stars, con- 
trary to expectati ons based on the 'dual outflow' sgB[e] model 
(lAret etal.ll2012h . Therefore, it is clear we have yet to arrive 
at a complete understanding of these objects. To systematically 
study the circumstellar environments of sgB[e] stars, we have 
conducted a CRIRES survey of Galactic sgB[e] and B[e] stars. 
The CRIRES observations presented here were drawn from this 
survey and demonstrate the quality and potential of the data. 
Analysis of the entire sample will be presented in a forthcoming 
paper (Muratore et al. in prep.) and follow-up VLTI observations 
are planned. 
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Fig. 5. A comparison between the velocity structure (upper row) and brightness distribution (lower row) of the best fitting disc (left 
column) and equatorial outflow (right column) models. The contours in disc velocity image are drawn at: 15, 30, 45 & 60 kms -1 
and the contour levels of the outflow velocity image are: 10, 30, 50 & 70 km s -1 . The maximum velocity shown in the outflow image 
is 100 kms -1 , but the velocity at the inner edge of the structure is ~400kms -1 . The brightness distributions are normalised to the 
contribution of the brightest element. 



6. Conclusion 

This paper presents high spectral resolution (25 & 6 kms -1 ) 
AMBER and CRIRES observations of the supergiant B[e] star 
HD 327083. The observations spectrally resolve the CO band- 
head emission of this object for the first time and provide new 
insights into its immediate environment. Here we list the key re- 
sults. 

- We find no evidence that the CO bandhead emission of 
HD 327083 is variable. Since the observations span the ma- 
jority of the estimated period of the binary system, this sug- 
gests that the structure of the CO emitting region is not 
strongly dependent on orbital phase. 

- The CRIRES observations of HD 327083's CO 1 st overtone 
bandhead emission are fit to an excellent degree by a model 
of a Keplerian rotating disc. 

- A differential phase signature corresponding to a photo- 
centre displacement of ~0.15 mas is observed. This is well 
reproduced by the same model that fits the high spectral reso- 



lution CRIRES spectrum. The agreement between the model 
and observations in both the spatial and spectral domains is 
compelling evidence that the CO emission of HD 327083 
originates in a Keplerian disc. 

- The inner edge of the best fitting CO disc is approximately 
3+0.3 AU. This places it interior to the previously detected 
continuum emitting disc (which has an inner radius of ap- 
proximately 5 AU), but is sufficiently large for the CO emit- 
ting disc to be circumbinary. 

- We show that a model of an equatorial outflow cannot re- 
produce the data as well as the Keplerian disc model. This 
suggests that the standard sgB[e] scenario is not applicable 
to HD 327083, which supports the hypothesis that the B[e] 
behaviour of HD 327083 is due to binarity. 
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